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Abstract 

The London ground-state energy formula as a function of number density for a system of identical boson hard spheres, 
corrected for the reduced mass of a pair of particles in a "sphere-of-infiuence" picture, and generalized to fermion hard- 
sphere systems with two and four intrinsic degrees of freedom, has a double-pole at the ultimate regular (or periodic, 
e.g., face-centered-cubic) close-packing density usually associated with a crystalline branch. Improved fluid branches are 
contructed based upon exact, field-theoretic perturbation-theory low-density expansions for many-boson and many-fermion 
systems, extrapolated to intermediate densities via Pade and other approximants, but whose ultimate density is irregular 
or random closest close-packing as suggested in studies of a classical system of hard spheres. Results show substantially 
improved agreement with the best available Green-function Monte Carlo and diffusion Monte Carlo simulations for bosons, 
as well as with ladder, variational Fermi hypernetted chain, and so-called L-expansion data for two-component fermions. 

PACS: 67.40. Db; 64.10.+h; 05.30. Fk; 05.30. Jp; 05.70.Ce 
Key words: boson and fermion hard-sphere fluids; random close-packing; liquid 4 He and 3 He; nuclear and neutron matter. 



1 Introduction 

An analytical formula for the ground-state energy E of an ./V-hard-sphere-boson system of volume f2 for all particle-number 
densities p = N/il was proposed by London [JJ as 

2Trh 2 c 11 

E N = TR T7T- 1 

771 (p~ X/3 -Po f ) 2 (p- 1/3 + bPo 1/3 ) 

where m is the particle mass, c is the hard-sphere diameter and the constant b equals 2 5 / 2 /tt — 1. Here, p = V^/c 3 is the 
assumed ultimate regular (or periodic) close-packing density at which a system of identical classical hard spheres close-pack 
in a primitive-hexagonal arrangement, e.g., face-centered-cubic or hexagonal. As remarked by Rogers [2] this is what "many 
mathematicians believe and all physicists know" to be the case. However, the Kepler 1611 conjecture [3] that p = \f 7 2.j(? is 
the ultimate packing density for identical hard spheres seems to be approaching theorem status [4] after many attempts of 
proof. 
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The justification given for ([T]) is that it reduces smoothly to limiting expressions at both low and high densities, namely 

E/N ► {2nh 2 /m)pc (2) 

p^O 

E/N > A(h 2 /2m)(p-^-p- 1/3 )- 2 (3) 

P^Po 

but gives no indication of a "freezing" or Kirkwood [5] phase transition at some number density p between and p . Here 
A = 7r 2 /2 1 / 3 ~ 7.8335 is a constant called the residue of the second-order (or double) pole at close packing. Using the 
polyhedron cell method suggested in Ref . [6] , the value of A has been predicted [7] theoretically to lie within the rigorous 
range 

1.63 <A< 27.0. (4) 

The low-density leading term j2]) is the celebrated Lenz [8] term, calculated by him as the leading correction to the energy 
arising from an "excluded volume" effect. The Lenz term has finally been rigorously established [9 . The limit Q comes 
from the lowest Schrodinger equation eigenvalue of a particle in a spherical cavity, and is just the kinetic energy of a point 
particle of mass m inside the cavity of radius r — c, where r is the average separation between two neighboring hard spheres 
and r = {\/2/ p) 1 ^ by assuming a primitive- hexagonal packing arrangement for the cavities. 

More recently it was found [TU], however, that the arguments leading to the high-density limit of the original [T] (boson) 
London equation (fTJ) are flawed by a fundamental error: the spherical cavity of radius r — c alluded to above in reality refers 
to the "sphere of influence" of two particles. Thus, the particle mass used in obtaining Q should refer to the reduced mass 
m/2 of the pair. This yields the constant 

fe = 2 3 / 2 /vr-l (5) 

instead of the constant 2 5 ^ 2 /tt — 1 given by London for (fl]). The result (fl]) with ([5]) is designated the modified London 
(ML) equation. It continues to satisfy © as this is independent of the constant b but the residue A in §5§ now becomes 
2 2 / 3 7r 2 ~ 15.667 instead of the previous 7r 2 /2 1 / 3 ~ 7.8335 associated with the original London equation, and fully agrees 
with the empirical residue of 15.7 ± 0.6, extracted by Cole [TT] from high-pressure crystalline- branch data in 3 He, 4 He, H2 
and D2 systems. Moreover, this ML equation exhibits dramatically better agreement than the original London (L) equation 
with Green-function Monte Carlo (GFMC) [T2] computer-simulation data points for both fluid and crystalline branches of 
the boson hard-sphere system. 

A generalized London equation has also been proposed [13] for N-fermion hard-sphere systems with v intrinsic degrees 
of freedom for each fermion. Here v = 2 for, say, liquid 3 He or neutron matter, both constituent fermions of which have spin 
1/2, and v = 4 for nuclear matter consisting of both neutrons and protons of spin 1/2. As v is essentially the maximum 
occupation in a given single-particle quantum state, it can be taken as infinite in the case of bosons. For fermions, two 
differences appear with respect to the boson London Formula: (a) unlike the boson case, the ground-state kinetic energy 
for fermions is nonzero and is added as a well-known [14j ^-dependent leading term; and (b) the constant b is allowed to be 
^-dependent, being replaced by 

6„(i/) = [(i/-1)/i/](6+1)-1 (6) 

which clearly approaches b as v — > 00. The latter form also ensures a v- independent energy at close-packing where, since 
the spheres can be labeled so that indistinguishability as well as particle statistics disappears, as expected in this classical 
limit. Substitution of b v for the constant b in (fTJ) gives a generalized form of the modified London equation (ML„) 

l\2nh 2 c 1 1 

E/1M = U v p-< " + I — 

ith 



E/N = c v P 2/3 + ( - — i\=^li : — : (7) 



3h 2 /6tt^ 2/3 



lUm \ v J v^oc 

For v — * 00 , b(y) — * b according to (6), and ([7]) goes over into the boson case (TTJ) because C v vanishes in this limit. The 
low-density limit of (0 is 

E/N^ ClJ p 2 ' 3 +( V —±Y-^pc (9) 
p->0 \ v ) m 
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where the second term on the rhs is the Lenz term for ^-component fermions in 3D. On the other hand, for p — > p Q = \[2/c 3 
one sees that ([7]) reduces to ([3]) as it should. In other words, hard-sphere fermions, bosons or "boltzons" must all close- 
pack regularly at the same density. From this it follows that the residue for bosons or fermions is the same and equal to 
2 2 / 3 7r 2 ~ 15.667, in excellent agreement with the empirical Ref. [TT] value of 15.7 ± 0.6. 

For bosons, in addition to the Lenz term for the low-density fluid branch, several higher-order corrections to the 
ground-state energy per particle have been derived using quantum field-theoretic many-boson perturbation theory [T5]|16|. 
They give 

E/N = ^^{l + C^pc 3 ) 1 ' 2 + C 2P c 3 ln(pc 3 ) + C 3 pc 3 + (pc 3 )} (10) 

for pc 3 <C 1, where C\ = 128/15y^r and Ci = 8(47r/3 — V%), but C3 is an as yet unknown constant. Here, c denotes the 
S-wave scattering length for a general potential; for a hard-core potential it is just the hard-sphere diameter. The series is 
clearly not a pure power series expansion, and is at best an asymptotic series. 
Similarly, for an ./V-fermion hard-sphere system the corresponding series is |17j 

E/N = ^^L{l + Cl (k F c)+C 2 (k F c) 2 

+ [C 3 r a /2c + C i A 1 {0)/c 3 + C 5 ](k F c) 3 + C 6 (k F c) 4 \n(k F c) 

+ [C 7 r /2 C + C 8 <(0)/c 3 + C 9 ](fc FC ) 4 + O (^c) 4 | (11) 

for k F c <C 1 and where the Cj (j = 1,2, ...,9) are dimensionless coefficients depending on v\ they are given in Ref. [T5] for 
v = 2 and v = 4. The Fermi momentum Tik F is defined through the fcrmion- number density 

p = N/n = vk 3 F /Qn 2 (12) 

with fl the system volume, so that the Lenz term expressed in terms of p is identical to the boson Lenz term apart from a 
factor of [y — \ )/v which is the average number of fermions the Pauli Principle allows a given fermion to interact with at the 
shortest possible range. 

Unfortunately, both low-density expansions (fTTJf and (fTTj) lack accuracy at moderate to high densities, including the 
saturation (or equilibrium, zero-pressure) densities of liquid 4 He (y — 00) [19] and liquid 3 He {v = 2) or nuclear matter 
(y = 4). However, one can extrapolate the series for hard-sphere systems to physical and even to close-packing densities 
through the use of Pade [20] and/or a modest extension of these called the "tailing" [21] approximants. The so-called quantum 
thermodynamic (or van der Waals) perturbation theory (QTPT) [22, 23J has provided fairly accurate representations of the 
fluid branch of the equation of state of quantum hard-sphere systems [24] . even beyond freezing (or, Kirkwood) phase 
transition densities, but without sufficient credibility as one approaches close packing. This is clear since one does not 
possess a single ground-state energy function with implicit information of both fluid and crystalline branches, with presumably 
different close-packing ultimate densities. 

In Section 2 we discuss the double- (or second-order-) pole behavior for the equation-of-state fluid branch conceivably 
ending at random closest close-packing, instead of the regular close-packing at which the crystalline branch terminates; in 
Sees. 3 and 4 we construct analytical expressions for the fluid branches for hard-sphere bosons and fermions, respectively. 
Sec. 5 gives our conclusions. 



2 Double-pole conditions at close-packing 

We shall assume that the fluid branch of the hard-sphere equation of state will terminate not at the regular close-packing 
density p but rather at the random closest close packing (rccp), sometimes called the Bernal, density p rccp (or "maximally 
random jammed" packing 25j). Its value was originally determined empirically 26J with actual ball-bearing packings. Near 
the density p rccp we expect, based on ([3]), that the energy for a hard sphere boson or fermion gas has the following behavior 

E/N ► A(h 2 /2m){p- 1 / 3 - p- c % 3 )- 2 (13) 
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with A the residue which could be different for each system. Random close-packing densities range [27) from about 0.06,Oq 
to 0.86p = p rccp . 

The derivative of (fT5|) with respect to p then tends asymptotically to 

d(E/N) A(2/3)(H 2 /2m) 

reap 

while 



dp P^ Prccp ( p -l/3 _ ^1/3)3^4/3 

d1n(E/N) 2/3 



(15) 



^ P-^W (p-1/3 _ p-V*)p*/Z 

is residue independent. We shall assume that ^4 is the same for boson as for fermion hard spheres and that their rccp density 
is likewise identical since at closest close-packing the particles become localized by definition, enabling one to formally 
label each particle; this makes them distinguishable thus rendering (quantum) statistics irrelevant. Note that the pressure 
P = p 2 [d(E/N)/dp] from (fT4]) also diverges as p — ► p rccp , as expected. 

3 Boson hard-sphere fluid 

In order to extrapolate the low-density series (fTT)|) to higher densities we start by writing it as 

E/N =^-pce (x) (16) 
m 

where x = (pc 3 ) 1 / 2 and 

e (x) = 1 + dx + C 2 x 2 lnx 2 + C 3 x 2 + 0{x 3 lnx 2 ) (17) 
for Alternatively, one can rewrite this series as 

e„ 1/2 (x) = l + K lX + K 2 x 2 lnx 2 +K 3 x 2 + 0(x 3 lux 2 ) (18) 

where the K^s are expressible in terms of the Cj's. As C 3 is to date unknown, consequently K 3 is also unknown. Values of 

^ 1/2 

the CiS and K^s are given in Table [T] We analyze the series e (x) instead of the series eo(x) to ensure that any zeros in 

I/O 

its extrapolants, say e (x), are double (or second-order) poles in the energy as one expects at any kind of close packing. 
The extrapolants are generated as a quotient of two polynomials such that on expansion one recovers the first terms of the 
original series. Series (fT5)) with three terms beyond unity has twelve extrapolants correctly generated in Ref. but fitted 
there to erroneous (i.e., to one-half the correct values) GFMC data points [12]. Adjusting various extrapolants [24j to best-fit 
the four known GFMC data points ensures a good value for the unknown coefficient K3 in (|18jl . The extrapolant labeled 
"XI (bosons)" in Fig. 2 of Ref. [24] had the least mean-square deviation with respect to the four GFMC fluid-branch data 
points. Therefore, we adopt it as our best initial extrapolant. The ground-state energy per particle for boson hard spheres 
was thus represented (symbol =) by 

E/N = —pce (p) (19) 
m 

with K3 ~ —27.956. However, as diffusion Monte Carlo (DMC) calculations became available [33] spanning a wider range of 
densities in the fluid region than GFMC data, we realized that although our expression XI(x) in Eq. (17) of Ref. [24] agrees 



Bosons (y = 00) 


i =1 


2 


3 


a 


4.81441778 


19.65391518 


"73.296" 


Ki 


-2.40720889 


-9.826957589 


"-27.956" 



Table 1: Coefficients C\ and Ki for bosons appearing in l|17p and (|18p. respectively. Numbers in quotation marks are 
determined as indicated in text. 
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well with DMC and GFMC data around the freezing transition, its disagreement with the DMC data at low to intermediate 

— 1/2 

densities suggested the possibility of improving the extrapolant. As will be seen, the new extrapolant e predicts a random 
closet close-packing (rccp) density p rccp / p — 0.776 which is only about 10% below the classical hard-spheres empirical [55] 
rccp value ~ 0.86 mentioned before and also assumed to be the ultimate rccp density for quantum hard-sphere fluids. 

In order to improve the fluid-branch expression of Ref. 124] for low to intermediate densities we use the two double-pole 
conditions (fl3|) and (fT4|) which lead to the following conditions on the extrapolant eo(x) to be used in (fl9|) . namely 

: 1/2 (p) ► [A/l-Kpc]- 1 ' 2 (p- 1 / 3 - p-^) ► 0. (20) 

P->Prccp P~*Pr-aap 



This is equivalent to 
The condition lfl4|) gives 



Clp P^Praap 



(21) 



Strictly, any log term should be accompanied by a constant, if known, because the scaling of p by c 3 is arbitrary. We 
thus propose the representation of eo (x) in (fl~7|) as given by 



-1/2^ j_ 1 + K 1 x + fix 2 + 7X 3 = ^_i/2 
1 — Ki x 2 In x 2 + ax 2 



(^ ^S^ ^^) (22) 



where a, (3 and 7 are to be determined from ([20]) and ^2TJ and by fitting both DMC [29] and GFMC (,12], Table I) data. In 
this approximant the terms in x 2 In x and x 2 are kept together. Condition ([20]) applied to ([22]) gives 



1 + Ki x rccp + (3x 2 rccp + 7x 3 ccp = . (23) 



The second condition ([2T]) can be rewritten as 



d(e B 2 ) dx 



K\ + 2(3x rccp + 37X 2 



dx dp 

Substituting (|23[) in the last equation we obtain j3 in terms of a, namely 

2 



TCC \ n = "(1/3) [A/^c]- 1 ' 2 p~X 6 ■ (24) 



1 K2 x 2 ccp In x 2 ccp + ax 2 ccp 2x rccp 



Px 2 rccp = 3 + 2K lXrccp - [A/An]^ 2 1^(1 - K 2 x 2 ccp In x 2 ccp + «x 2 ccp ) . (25) 



Now substituting ([25]) in ([23]) we arrive at 

2 

rccp — " t "l-Tccp — 1/1 j g^rccpV^ — "2 -^ccp lli '''rccp t "-^rccp 

Introducing ([25]) and ([26]) in ([22]). we get 



l X rccp — ^ + "^1 X r ccp [A/47r] ^ n X rccn(^ ^2 X rccn^ n X rccn "F aX rccn) ■ (^6) 



£ 0B 2 ( X ) 



1 + Jfi x + (x/x rccp ) 2 {-3 - 2^!X rccp } + (x/x rccp ) 3 {2 + ifix rccp } + L4/4vr] 1/2 -x^ 3 p 

x (1 - K 2 x 2 ccp In x 2 ccp + ax 2 ccp ){(x/x rccp ) 2 - (x/x rccp ) 3 }] [1 - K 2 x 2 In x 2 + ai 2 ]' 1 (27) 

from which after some algebra one obtains a single equation for a, namely 

e os /2 ( a '^' x ) (! _ ^2X 2 In x 2 ) - 1 - K\ x + (x/x rccp ) 2 {3 + 2K x x rccp \ 
2 

[-i4/47r] g^rccpl^ ^-2 %rccp ^ **Vccp){ (*^/ ^Vccp) / %rccp) J / %rccp) ~h i^C^X rcC p} 

= a{x 2 ccp [A/47r]- 1/2 -x r £? p [(x/x rccp ) 2 - (x/x rccp ) 3 ] - x 2 [e^ /2 (a, A,x)]} . (28) 

— 1/2 

where we have explicitly written the dependence of e 0B (x) on a and A. To determine a from the DMC [29] and/or GFMC 
data we must calculate the values af MC [from ([28]) after replacing eQg 2 (a, A, x) by the ^q1dmc( x F MC ) obtained from 
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(PI) as {2Trh 2 pcN/Em) 1 / 2 with £/iV the energy from DMC calculations] for each xf MC for i = 1, 2, • • • N values, and then 
minimizes Y2i=ii a F MC — ct) 2 by imposing 



, n 



which gives 

N 



i=l 

Since the fluid branch GFMC data are a subset of DMC data, we have used these to calculate a here, determining A in the 
next step. For residue A fixed at 2 2 / 3 7r 2 ~ 15.667 as described below (0, we obtain an optimal a ~ 114.282 which from ([25]) 
and (gnj) leads to (3 _ 74.0891 and 7 ~ -65.9475. The curve then corresponding to is labeled Bl in Fig. [TJ 

Alternatively, if we allow the residue A to be free one may ask for a solution minimizing Eili [ £ o-dmc (■'■f — 
e g 2 (a, A, xf MC )] 2 with respect a and A, i.e., 

A N N A 

^ Y>*-nMC^ MC ) ~ ^(a, ^ x? MC )f =-J2 ^Jomc^ MC ) - A, z?»°)]±e£ /2 (a, A, X ™°) = 

i—l i—1 



or 

AT 



I. 2 , /j.l/r vl, Xi)x 2 + (A/47r) 1 / 2 |x r ccp3 : rccp[( ;z '/ :z; i , ccp) 2 (^/^rccp) 3 ] 



^o-DMC^i ) £ os ^,*< jj ________ _U ^yj 

as well as of 

, iV iV , 

lVtr 1/2 ( T DMC\ _ .-VS/ , r i)MC\l2 _ _V^ 2 r -1/2 /DMC\ -1/2 / , DMC^i a -1/2 , , DMC\ — q 

i=l i=l 

or 

Er,r -1/2 / DMC\ -1/2/ a DMC\[ (^V^ 71 ") 3 ^ 2 J^Xrccpjl — -K2 ^rcep m ^rccp + aa: rccp) [(^/^rccp) 2 ~ (^/Zrccp) 3 ] _ 
1 — K2 x 2 In x z + ar 

i=l 

(30) 

Under the two conditions (|29|) and (|30|) we find an optimal A ~ 11.8715 and an optimal a ~ 169.516, leading to (3 ~ 124.1 
and 7 _ —111.296. This procedure gives the curve labeled B2 in Fig. Q] Note that the residue 11.8715 is now being associated 
with the random closest close-packing (rccp) density 0.86/o of hard spheres. This value of A is somewhat smaller than the 
residue 15.667 at regular close-packing density p , though still within the rigorous range stated in gl). Figure [2] compares 
the previous fluid branch expression XI(a;), Eq. (17) of Ref. [53], with the present extrapolant (f2"2"]l labeled B2, both as full 
curves. The dashed curve is the modified London (ML) formula (Q]) that connects smoothly with the crystalline branch. 
Open circles and squares are GFMC data for fluid and crystalline branches, respectively. Dots represent DMC [29 data 
spanning a wider range of densities in the fluid region than the GFMC data. The new expression B2 shows dramatically 
better agreement with DMC data for intermediate densities, as well as agreeing well with both DMC and GFMC data around 
the freezing transition mentioned in Table I of Ref. [12 . Fig. 3 is an enlargement of Fig. 2 at low densities to show the 
remarkable agreement of B2 with the DMC data. 

4 Fermion hard-sphere fluid branch 

The ground-state energy per particle for fermion hard-sphere fluids (fTT]) can be written as 

E/N = llJ&eoix), x = k F c (31) 
5 2m 
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0.0 0.1 0.2 0.3 0.4 

P/Po 



Figure 1: The quantity e ^ 2 = 2-kTi 2 p c N / m E = [1 — (p/ Po) 1 ^ 3 ]\/l + b (p/Pa) 1 ^ 3 as a function of x/xq for boson hard 
sphere systems: Bl and B2 refer to (j2"2"]l and (|2~4"]) with A ~ 15.7 and A ~ 11.9, respectively. Larger dots are GFMC fluid 
data and smaller dots refer to DMC (fluid) calculations. 

with 

e (x) =1 + d x + C 2 x 2 + (C 3 /3 + C4/3 + C 5 ) x 3 

+ C 6 x i lnx + (C 7 /3-C 8 /3 + C 9 )x i + o(x i ) (32) 

for x = kpc I, p = N/Cl — iskp/6ir 2 being the number of fermions N in the enclosed volume fl. We shall examine both 
v = 2 (corresponding to liquid 3 He and neutron matter) and v = 4 (corresponding to nuclear matter). 

4.1 Fermions with v = 2 

For zv = 2, C6 = [17] so that (|32|) simplifies to the pure power series 

e (x) = 1 + D 1 x + D 2 x 2 + D 3 x 3 + D 4 x 4 + o(x A ) (33) 

where the C^'s have been determined in terms of the ZVs. As in the boson case, instead of eo(x) we consider the series 

e~ 1/2 {x) = 1 + Fi x + F 2 x 2 + F 3 x 3 + F 4 x 4 + F 5 x 5 + o{x 5 ) (34) 

where the Fi's depend algebraically on the ZVs in a simple manner, F$ being unknown. Values of Di and Fi are given in 
Table [2j We use this simple power series to construct the usual Pade extrapolants. The approximants to (f34l) with four 
terms beyond the trivial unity were analyzed in Ref. [30] where it was concluded that the best approximant was the Pade 
[0/4] (x). However, this function does not have a zero in the region of physical interest, i.e., < p/p < 1, which implies 
< x = k F c < 3.47 since p — kp/3n 2 . Accordingly, the energy does not manifest a close-packing density as it should. This 
deficiency made it advisable to introduce the fifth term F 5 x 5 in (|34[) . Although in Fig. 1 of Ref. [3T] only five of the six 
two-point Pade approximants [L / / M](x) with L + M = 5, are shown, all six approximants were analyzed here to adjust F5 
so as to ensure a zero associated with a random close-packing in the physical region. The approximant eo (x) and the position 
of its zero were chosen in such way that the QTPT applied in Ref. [30] to calculate the ground-state energy of 3 He with 
the Aziz interatomic potential [32] reproduces the corresponding GFMC [34] data. (In this treatment, the Aziz potential 
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1.0 




P/Po 

Figure 2: The quantity e = \J 2-rtTi 2 p c N / m E = [1 — (p/ Pq) 1 ^ 3 ] y/l + b (p/ Po) 1 ^ 3 as a function of pj p a for the boson hard 
sphere system: XI is the fluid branch approximant of Ref. [53], Fig. 2; B2 refers to (f2"2")l and (|2~4"]) with A ~ 11.9; ML is the 
modified London formula ((T|). Open circles and squares are GFMC data for the fluid and crystalline branches, respectively, 
and dots are DMC data points. 



was decomposed via the well-known Barker- Henderson (BH) [33] scheme as described in Ref. [30].) Eventually, the best 
extrapolant was found to be the two-point Pade approximant 

\ • ta/fov ^ _ N Q + N x x + N 2 x 2 + N 3 x 3 _ -1/2, v 
e n (x)= 3//2 (#) = — — — — t. = e n (x) (35) 

where 

AT = F 2 F 4 -F* 

Ni — F^F^ + F1F2) — F2F5 — FiF 3 

N 2 = (F 3 - F 1 F 2 )F 5 - Fi + (FiFs + Fi)F 4 - F 2 Fi 

N 3 = (F 1 F 3 -F 2 2 )F 5 -F 1 Fi + 2F 2 F 3 F i -F$ 

Mq — F2F4— F 3 

Mi — F 3 F 4 — F 2 F 5 

M 2 = F 3 F 5 ~F 2 . 

The extrapolant (|35| satisfies [3/ /2](x = 3.13) = 0. Hence the ground-state energy per fermion for v = 2 becomes 

E/N = |^i{[3//2](x)}- 2 (36) 



v = 2 


i =1 


2 


3 


4 


5 


A 


0.353678 


0.185537 


0.384145 


-0.024700 


"-0. 265544" 


Fi 


-0.176833 


-0.045863 


-0.156677 


0.109672 


"0.130830" 



Table 2: Coefficients Di and Fi for ^ = 2 appearing in (|3"3"|) and (j3"4"]) , respectively. Numbers in quotation marks were 
determined as indicated in text. 
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71 1 1 

0.0001 0.001 0.01 0.1 

P/Po 



Figure 3: Enlargement of Fig. Oat low densities. 



with a random closest close-packing density p rccp / p = 0.732 only 15 % smaller than the empirical [SB] value p rccp /p — 0.86. 
The coefficient F& is listed in Table O in quotation marks. In Fig. 2] we show the expression 

e - 1/2 = [3h 2 (6Tr 2 p/^ 2 / 3 N/l0mE}^ 2 = 1 + 

[20^ - l)/3H(2 1 / 4 ^/6^) 2 / 3 {[(p/p )- 1 / 3 - l] 2 [(p/p )- 1/3 - bWip/pof/ 3 }- 1 (37) 

as a function of p/p for fermion hard spheres. Here b(y) is as defined in ©. For = 2 the fluid branch [3/ /2] (full curve) 
given by (|35| is close to the Ladder [35] (open squares), the variational Fermi hypernetted chain (VFHNC) [36] (plus-sign 
marks) , and the so-called L-expansion data [37] [38] (open triangles) . Fig. [4] shows good agreement over the entire range of 
available data. 




Figure 4: The expression ([37]) as function of p/p for fermion hard spheres with v = 2 labeled [3//2] and with v — 4 labeled 
XII (full curves). Dashed curves are the corresponding modified London ML^ formulae, but note that the ML2 dashed curve 
almost coincides with the full curve [3//2]. 
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In order to improve the many-fermion ground-state energy equation of state we include the next term in 

I/O 

which is then used to generate all Pade approximants of order six to the series e ' (x). The lack of a logarithmic term 
x A In a; is due to the Pauli principle [17] . Such a term arises when there are three independent hole lines. But for v = 2 there 
can be at most two lines of the same spin. Thus the Pauli principle reduces the size of the term by a factor of the density. 
We thus expect the first such term for v = 2 to be 0(x In a;). The unknown coefficients Fq and F$ are determined from the 
two double-pole conditions (fl~3[) and (fT4|) . which become 



e Q 1/2 (x) ► {l-x/x rccp )[hA/Z(?,K 2 f' z ]- 1/2 (38) 

and 



eo{x) + *^l > 5A/3(3tt 2 ) 2 / 3 

2 dx p^p rccp {l~x/x rcC p) 3 



ith A ~ 15.667. 




Figure 5: Improved extrapolants for the many-fermion hard-sphere gas with v = 2. 



For each Pade approximant of order six we determined F 6 and F 5 as shown in Table Approximants [4/2] and [0/6] 
did not exhibit the double-pole conditions. The other four approximants are plotted in Fig. [5] together with the Ladder [55] 
(open squares), the variational Fermi hypernetted chain (VFHNC) [36) (plus-sign marks), and the L-expansion [37j (open 
triangles) data for v = 2, from which we conclude that the approximant [3/3] (x) is the best. Fig. [6]is a semi-log enlargement 
of Fig. [5j In Fig. Owe compare both the new improved expression [3/3] (x) and the previous best energy expression, i.e, the 
two-point Pade approximant [3//2](a;) reported in Ref. [53] and supported by Ladder, VFHNC and L-expansion data. 
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Pade 


F5 


F6 


[5/1] 


-0.0272548 


0.0038205 


[4/2] 


-.20 


no solution 


[3/3J 


-0.0130625 


U. 0039120 


[2/4] 


-0.0395076 


0.0415222 


[1/5] 


-0.0115902 


0.01887153 


[0/6] 


-0.1276 


no solution 



Table 3: The F5 and Fq coefficients for v = 2 that follow from conditions (|38|) and (|39| for all sixth-order Pade approximants 
with residue A ~ 15.667 and random closest close-packing density /0 rccp = 0.86p . 




Figure 6: Enlargement of Fig. [5] at low densities. 

4.2 Fermions with i> = 4 

For fermions with v = A p2|) becomes 

e (x) = 1 + £>i x + D 2 x 2 + D 3 x 3 + D 4 x A In a; + D 5 a; 4 + o(a; 4 ) (40) 
for x = fci?c <C 1 and we recall that p = vkp/Gir 2 . As for bosons or for fermions with v = 2, we analyze 

eo " 1/2 (a;) = l + Fi x + F 2 x 2 + F 3 x 3 + F 4 x 4 \nx + F 5 x 4 + o(x 4 ) (41) 

with all Fi (i — 1,2,3,4) known. Values of Di and Fi are given in Table 2] Unlike the v — 2 case, this series is not a pure 
power series as it contains logarithmic terms. Its so-called "tailing" [5T] approximants are given in Table III of Ref. 28 . 
Of all the possible approximants using only the known coefficients, only the forms II and XII are free from flaws and have 
residues within the bounds ([4]). Of these two forms, II has a residue less than that predicted in Ref. [7]. Hence we chose 
form XII, which is plotted in Fig. 0] as the full curve labeled XII. 

In this case E/N can be written as 

3h 2 k 2 



E/N=—^e (x) (42) 
10m 



where the series (|4ip is represented as 



-1/2, . . YTT/ . 1 + (Fx - F 3 /F 2 )x + {F 2 - F 1 F 3 /F 2 )x 2 _ 1/2 
so (x)=XH(x) S = £q ( 
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Figure 7: Comparison of quantity ([37]) as function of p/p for many-fermion hard spheres with v = 2, for the the previously 
best approximant [3//2](x) |24j and the new improved one [3/3] (x), full curves. Dashed curve is modified London formula. 



v = 4 


i =1 


2 


3 


4 


A 


1.061033 


0.556610 


1.300620 


-1.408598 


Ft 


-0.530517 


0.143867 


-0.5806558 


-0.704299 



Table 4: Coefficients Di and -Fi for v = 4 appearing in (1401) an d (|4Tj) . respectively. 

We also plot the corresponding VFHNC data (plus-sign marks) and L-expansion data (open triangles). In terms of energy, 
our results are slightly below the VFHNC points, with agreement improving at lower densities. On the other hand, the XII 
approximant lies just above the L-expansion data over the range of densities where data are available. 

In order to improve the v = 4 many-fermion hard-sphere ground-state energy equation of state, the energy series ()llj) 
was written as 

I - r£f - - l^ D < + ^ + D ^ + D ** 3 + D - 3 

+L> 6 x 4 lnx + L> 7 x 4 + ---] (43) 
where x = kpc and p = vkp/Gn 2 . The suggested representation for eo(x) here is 

e (x) = Di + D 2 x + D 3 x 2 + D 4 x 3 In x + D 5 x 3 + D 6 x 4 In x + L> 7 x 4 H (44) 

which leads to 

e (x)~ 1/2 = F\+ F 2 x + F 3 x 2 + F 4 x 3 \nx + F 5 x 3 + J F 6 x 4 lnx + F 7 x 4 + • • • (45) 

with D\ to D4 known and equal to the values given in the Table 2J The coefficients F\ to F4 are different from those in 
Table SI but they are derived simply from the ZVs and so are also known. They are 

F 1 = 1/VA; F 2 = -D 2 /2D\ /2 - F 3 = {3D 2 2 - AD^) I%d\' 2 ■ F 4 - -8£>?£> 4 /16£>[ /2 ; 
Fs = (-5-Dl + 12DxD 2 D 3 - 8D 2 D 5 )/16Dl /2 ; F 6 = 32L> 2 (3L> 2j D4 - 2D 1 D 6 )/128D 9 1 /2 . 
We have also investigated the representation 

-1/2, , F 1 + F 2 x + F 3 x 2 + bx 3 



1- (F 4 /Fi)x 3 lnx 
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for which the two double-pole conditions (JT3J) and ([Hp imply that 



rccp 



[Fi 



F 3 x 2 



rccp} 



and 



2F 3 x r 



3bxt 



5.i" 



rccp 



3 

717/ 



1 



£4 



*^rccp lri *^ r ccp ^ CLCC 



3 

rccp 



The values of a and 6 so determined are —0.0924883 and 0.171942, respectively. This representation is unsatisfactory because 
it has what applied mathematicians call a "defect." Unfortunately it is in the physical region < x < x rccp . The problem is 
not uncommon and stems from a pole and a zero lying very close to each other. 



5 Conclusions 

Based on known terms of field-theoretic perturbative low-density expansions we have constructed closed-form analytical ex- 
pressions as functions of particle density using Pade and other approximants for the energy per particle of the fluid branches 
of both many-boson and many-fermion quantum hard-sphere systems. Improvements with respect to previous work (notably 
but not exclusively that of Ref. [24]) have been achieved by assuming i) that the classical random closest close-packing 
hard-sphere densities are the ultimate fluid densities at which the energy diverges with a second-order pole and ii) proposing 
and imposing a value for the residue at the pole that is the same for either bosons or fermions as closest close-packing is 
approached and the hard spheres become distinguisable. Implementing these two conditions and taking advantage of recent 
diffusion Monte Carlo simulation data has allowed us to incorporate an additional term in the low-density expansion beyond 
that employed in Ref. 24J. The resulting determination of the best approximants has produced decidedly improved results 
for bosons as well as for two-component fermions, but not for four-component fermions. 
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